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ABSTRACT 



Context. High levels of deuterium fractionation in gas-phase molecules are usually associated with cold regions, such as prestellar 
cores. Significant fractionation ratios are also observed in hot environments such as hot cores or hot corinos, where they are believed to 
be produced by the evaporation of the icy mantles surrounding dust grains, and thus are remnants of a previous cold (either gas-phase 
or grain surface) chemistry. The recent detection of DCN towards the Orion Bar, in a clump at a characteristic temperature of 70 K, has 
shown that high deuterium fractionation can also be detected in PDRs. The Orion Bar clumps thus appear as a good environment for 
the observational study of deuterium fractionation in luke-warm gas, allowing to validate chemistry models in a different temperature 
range, where dominating fractionation processes are predicted to be different than in cold gas (< 20 K). 

Aims. We aimed at studying observationally in detail the chemistry at work in the Orion Bar PDR, to understand if DCN is produced 
by ice mantle evaporation, or is the result of warm gas-phase chemistry, involving the CH2D+ precursor ion (which survives higher 
temperatures than the usual HoD^ precursor). 

Methods. Using the APEX and the IRAM 30 m telescopes, we targetted selected deuterated species towards two clumps in the Orion 
Bar. 

Results. We confirmed the detection of DCN and detected two new deuterated molecules (DCO* and HDCO) towards one clump in 
the Orion Bar PDR. Significant deuterium fractionations are found for HCN and H2CO, but a low fractionation in HCO^. We also 
give upper limits for other molecules relevant for the deuterium chemistry. 

Conclusions. We argue that grain evaporation in the clumps is rather unlikely to be dominant, and we find that the observed deuterium 
fractionation ratios are consistent with predictions of pure gas-phase chemistry models at warm temperatures (T ~ 50 K). We show 
evidence that warm deuterium chemistry driven by CH2D^ is at work in the clumps. 

Key words, interstellar medium - astrochemistry - deuterium chemistry - gas/grain interaction 



1. Introduction 

, Despi te the low de uterium abundance in the universe (D/H ~ 
■ 10" , iLins k^ l2003h . high abundances of deuterated molecules 
' have been observed in dark clouds and star forming re- 
gions in the last few years, with detections of molecules 
conta i ning up to three ato ms of deuteriu m (ND3: iLis et al.l 
I2OO2I: lvan derTaket"ani2002. and CD3OH: iParise et al.ll2004h . 
with noteworthy fractionation effects (CDaOH/CHaOH- 1%, 
iParise et al.ll2004l) . 

Formation of such highly-deuterated molecules requires spe- 
cific physical conditions, which makes them very interesting 
probes of the environments where they are observed. In molec- 
ular clouds, deuterium is mainly locked into molecular HD. 
Efficient transfer of deuterium from this reservoir at the low tem- 



perature of dark clouds has to occur through ion-molecule reac- 
tions, and it has long been known that only a few ions react fast 
enou gh with HP to c ompete aga inst electron recom bination: H3 , 



CH+ (lHuntressll977h and C2H+ (iHerbst et al.ll987h . The deuter- 
ated isotopologues of these three ions are thus considered to be 
the precursors of deuterium fractionation in the gas phase. The 
transfer of deuterium from the HD reservoir to other molecules 
is initiated through the following exothermic reactions: 



H- HD - 
CH+ + HD 



H2D+ + H2 
CH2D+ + H2 



(1) 
(2) 
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C2H^ -H HD ^ C2HD+ + H2 (3). 

H2D+, CH2D+ and C2HD+ then transfer their deuterium to the 
other species through ion-molecule reactions. Exothermicities 
are 232 K fc erlich et al. 2002) for reaction (1), ~390K 
Asvanv et alJM04) for reaction (2) and ~550K (IHerbst et all 



1 19871) for reaction (3), so that the reverse reactions are inhibited 
at low temperatures. Efficient transfer of deuterium to molecules 
has been widely observed in cold regions where high levels 
of CO depletion are present, such as dark clouds or prestellar 
cores. In these environments, H3 is predominantly responsible 
for the fractionation. Reaction (2) is thought to be predominant 
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at slightly higher temperature (T~ 30-50 K), when (1) is not ef- 
ficient anymore due to an increased importance of its reverse 
reaction. Although the chemistry involving H2D^ is now basi- 
cally understood, thank s to the numerous detect ions of HaD^ 
in prestellar cores (e.g. ICaselh et al.ll2003L 2008^ and the par- 
allel development of chemical models (e.g. Rob erts et al. 2003; 
IWalmslev et al.l2004l:lFiower et alj2004tlPagani et a l. 2009), the 
contribution of the CH2D^ chemistry has so far not been obser- 
vationally investigated, due to the lack of observations target- 
ting intermediate temperature sources, warm enough so that the 
CH2D^ contribution becomes significant relative to H2D^, and 
cold enough for the chemistry not to be dominated by ice evap- 
oration. 

Recently, during an unbiased spectral survey of the Orion 
Bar using the APEX telescope, we detected a deuterated 
molecule (DCN) for the first time in a m olecular clurnp in a 
Photon-Dominated Region (hereafter PDR, iLeurini et al ]l2006l) . 
This was however unexpected due to the high temperature 
(T ~ 70 K) characteristic of this clump. The fractionation ratio 
deduced from these observations is 0.7 %, a value intermediate 
between the one observed in dar k clouds (L134, 5% , Turner 
1200 ll) and hot cores (O rion, 0.1%.ISchilke et al.ll 19921: a sample 
of hot cores, 0. l-0.4%. lHatchell et al.lll998l) . 

DCO^ was not detected in this survey, with an upper limit 
on the DCOVHCO+ ratio of ~ 0.1 % (see below), although ob- 
servations in other environments point to similar fractionation 
ratios for the two species (dark cloud L134: Tine et al. 2000, 
Turner 2001; low-mass protostar IRAS 16293: van Dishoeck 
et al 1995). We proposed that this may be an indication that 
chemistry involving CH2D^ as precursor for deuterium trans- 
fer to molecules is at work in the Orion Bar, making it a refer- 
ence environment to study in further detail reactions involving 
these routes. This possibility was confirmed by the theoretical 
modelHng study of Roueff et al. ( 2007). Since the DCN detec- 
tion bv iLeurini et al. (2006), DCO^ was detected towards the 
Horsehead PDR fPety et aL ,20071 in a cold (10-20K) conden- 
sation shielded from the UV illumination. The observations we 
present in this paper are targetting warmer regions than the con- 
densation observed by Pety et al. (200^. 

We present here a detailed investigation on the deuterium 
chemistry at work in the dense clumps in the Orion Bar, based 
on observations with the APEX and IRAM 30m telescopes. 
Observations are described in section 2, the physical conditions 
in the clumps (temperature and H2 densities) are derived in sec- 
tion 3, the relative abundances and fractionations for targetted 
molecules are determined in section 4, and the chemistry is dis- 
cussed in section 5. 



2. Observations 

The Orion Bar, located in the Orion A molecular cloud, is a well 
studied PDR, mostly because of it s nearly edge -on morphology. 
Its distance is estimated to 414 pc dMenten et al., 20 07). It is illu- 
minated by the young trapezium stars, located some 2' to the 
north-west. Previous studies have shown that this PDR has a 
clumpy structure (see s ection [3]l. The interferometric observa- 
tions of iLis & Schilkd (l2003h displayed a series of molecular 
clumps, as traced by the H'^CN(l-O) transition, located in the 
cloud some 10" behind the H/H2 tra nsition. In the presen t study 
we use the nomenclature defined by lLis & Schilkd (l2003h to re- 
fer to the clumps. 

We present in this section the APEX and IRAM 30m obser- 
vations of the Orion Bar, targetting particularly clump 1 and 3 of 



iLis & Schilkd (l2003h . Unless otherwise stated, the transition fre- 
quencies of the molecules are taken from the Cologne Database 
for M olecular Spectroscopy (Mulleretal. 2001; Miilleretal] 
l2005h . Table[T]summarizes the observations presented in this pa- 
per. 

2. 1 . APEX observations 

Using the APEX telescope on Chajnantor (Chile), we mapped 
the Orion Bar in the DCN(4-3) and H'^CN(4-3) transitions. 
The double-sideband APEX2a receiver (R isacher et al. 20061) 
was tuned to 289.0000GHz (DCN) and 345.3397G Hz (H'^CN) 
and connected to the two units of the EFTS backend jKlein et al] 
2006), each with 8192 channels, leading to a velocity resolution 
of 0.13 and 0.11 km/s respectively over the two times 1 GHz 
bandwidths. The APEX beamsize is 21" (respectively 18") at 
289 GHz (respectively 345 GHz). 

The (0",0") position of the map is q'(2000)=05''35'"25.3', 
5(2000)=-05°24'34.0", co rrespon ding to the "Orion Bai- 
(HCN)" position of Schi Uce et al.l ( 2001). the most massive 
clump seen in H'^CN (Lis & Schilke'7003), as well as the target 
of the spectral survey of Leurini et al. (2006). The maps were ob- 
tained using the on-the-fly mode, with a dump every 6". The ref- 
erence position was taken at the (600", 0") offset position from 
the center of the map. 

The observations were performed between July 19* and 
August 2"^ 2006, under very good to good weather conditions 
(with a precipitable water vapor ranging from 0.3 to 1.5 mm). 
The typical DSB system temperatures were 115 and 200 K at 
289 and 345 GHz respectively. 

Several CH30H(6-5) transitions were also present in the 
DCN(4-3) setup. DCN(5^) and HNC(4-3) were observed to- 
wards the two clumps on June 28*, 2007, with Tjys around 
200 K. 

Observed intensities were converted to Tmb using 
Tmb - T*/77nib where r/mh - 0.73 (iGiisten et al.ll2006l) . 

We focus here on the observations towards clumps 1 and 3, 
and analyze the spatial distribution across the Bar in a later paper 
(Parise et al. in prep). 

2.2. IRAM 30m observations 

We targetted the Orion Bar during three observing runs at the 
IRAM 30 m telescope. During the first run we made use of the 
ABCD receivers, targetting only the two clumps in the Bar. The 
second run consisted of mapping the Bar with the HERA re- 
ceiver Some complementary data on the two clumps were also 
acquired as part of a third run. 

2.2.1. Single pixel receivers 

Using the IRAM 30 m telescope, we observed different species 
toward the two brightest H'^CN clumps of the Orion Bar — 
"clump 1" at offset position (0", 0") a nd "clump 3" at position 
(-50", -40"), as denoted byE&Schilke (2003). Besides ob- 
serving different transitions of DCN to constrain the excitation, 
and looking for a lower excitation line of DCO^ than the one 
that was not detected with APEX, we selected for our search 
molecules that can be synthesized in the gas phase via channels 
invol ving the CH2D^ ion. Such molecules include HDCO and 
C2D dTurneil bOOl). We also searched for CH2DOH and HDO, 
in order to constrain any possible ice chemistry contribution. 
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The observations were performed from September 29th to 
October 7*, 2006, under variable weather conditions. Four re- 
ceivers were used simuhaneously to observe two different fre- 
quency bands (either in the AB or CD setup). The observed Hnes 
are Hsted in Table[T] The receivers were connected to the VESPA 
correlator in parallel mode, leading to different velocity resolu- 
tions depending on the transition. Addi tional ly, the Bonn Fourier 
Transform Spectrometer ( Klein et alj |2006|) was connected to 
two of the receivers, providing a 850 MHz bandwidth. 

Some additional data (DCN(3-2) as well as integration on 
CH2DOH towards both clumps) were taken during a third ob- 
serving run in May 2008, under poor (CH2DOH data) to moder- 
ate (DCN data) weather conditions. 

2.2.2. HERA observations 

The Orion Bar was mapped in selected methanol and formalde- 
hyde transitions using the HERA receiver, a heterodyne array 
consisting of two arrays of 3x3 pixels with 24" spacing. The ob- 
servations were performed during the winter 2007 HER A pool 
observing period. The full dataset is presented in Leurini et al.l 
( I2OO9I) . Here we analyse the methanol observations towards the 
two clumps to derive their physical properties. 

All intensities of observations from the IRAM 30 m tele- 
scope were converted to Tn,b using Tmb - |^T* where Bgff is 
the main beam efficiency, and Fefj is the forward efficiency. The 
main beam efficiencies decrease from 78% to 50% between 87 
and 241 GH43. Forward efficiency is 95% at 3 mm, 93% at 2 mm 
and 91% at 1.3 mm. 



3. The clumpy morphology of the Orion Bar 

The Orion Bar was shown to have an heterogeneous struc- 
ture, with clumpy molecular cores embedded in an inter- 
clu mp gas. The two-c omponent morphology was firs t infeiTed 
by iHogerheiide et al.l (1995) and I Jan sen et al] (Il995h . because 
single-dish observations of CS, H2CO and HCO^ could not 
be fit with a single density component. The clumpy structure 
was later confirmed di rectly by interferometric observations 
(lYoung Owl et al.ll2000l) . The clumps have a density of several 
10^ c m~^ while the density of the interclump gas is ~ 10^-10^ 
cm"^ (lYoung Owl et al.y 2000). Interferometric maps have shown 
that H'^CN is mostly confined to the clumps (Lis & Schilk^ 
I2003h . which are relatively cold (~70K) compared to the inter- 
clump gas (of typical temperature of 150 K). 

In the following, we discuss the physical parameters of the 
clumps, based on new observations. We first derive the H2 col- 
umn density (Section [3Tl ), and then the temperature and H2 den- 
sity based on methanol observations (Section [J!2] |. 

3. 1. The H2 column density of the clumps 

In this section, we attempt to derive the H2 column density in 
the two clumps as accurately as possible, in order to be able to 
compute molecular abundances relative to H2 in the following 
sections. This will then allow us to compare the measured frac- 
tional abundances to predictions of chemical models. 

3.1.1. Clump 1 



Table 1. Summary of the observations. 



Transition 


Frequency 


Telescope 


Beamsize 


Targetted 




(GHz) 






Clump 


DCN(2-1) 


144.8280015 


IRAM 30 m 


17" 


1, 3 


DCN(3-2) 


217.2385378 


IRAM 30 m 


11" 


1, 3 


DCN(4-3) 


289.6449170 


APEX 


21" 


map 


DCN(5-4) 


362.0457535 


APEX 


17" 


1, 3 


DC0+(2-l) 


144.0772890 


IRAM 30 m 


17" 


1, 3 


HDCO(2h-1,o) 


134.2848300 


IRAM 30 m 


18" 


3 


CH2DOH(2~l) 


89.3 (band) 


IRAM 30 m 


28" 


1, 3 


HDO(2u-2i.2) 


241.5615500' 


IRAM 30 m 


10" 


3 


C2D(2-1) 


144.3 (band) 


IRAM 30 m 


17" 


3 


DNC(2-1) 


152.609774 


IRAM 30 m 


16" 


3 


H'^^CN(l-O) 


86.3399215 


IRAM 30 m 


28" 


1, 3 


H'3CN(3-2) 


259.0117978 


IRAM 30 m 


9.5" 


1, 3 


H'3CN(4~3) 


345.3397694 


APEX 


18" 


map 


H'3CO+(1~0) 


86.7542884 


IRAM 30 m 


28" 


1, 3 


HCO+(1-0) 


89.1884957 


IRAM 30 m 


28" 


1, 3 


Hi3cO(2„-l,o) 


146.6356717 


IRAM 30 m 


17" 


3 


C"O(l-0) 


112.36 


IRAM 30 m 


22" 


3 


C"0(2-l) 


224.71 


IRAM 30 m 


11" 


3 


CH30H(5-4)2 


241.8 (band) 


IRAM 30 m 


10" 


map 


CHiOH(6-5) 


290.1 (band) 


APEX 


21" 


1,3 


CH30H(li-lo A) 


303.367 


APEX 


21" 


1, 3 


HNC(4-3) 


362.63 


APEX 


17" 


1,3 



' IPL database. 

^ HERA observations (see Leurini et al. (subm) for more detail). 



' See |http://www.iram.es/IRAMES/telescope/telescopeSummary/| 
telescope_summary.html 



An H2 column density of 9 x 10^^ c m~^ averaged on the 18" 
beam was derived towards clump 1 bv iLeurini et al.l ([2006) from 
analysis of the C'^0(3-2) emission line, assuming a rotational 
temperature of 70 K. 

We can also get an estimate of the H2 column density in the 
clump from the analysis of the dust emission observed in the 
frame of a project targetting clump 1 with the Plateau de Bure 
interferometer in March, A pril and December 2004 (follow-up 
project of the work from lLis & Schilkell2003l) . The observations 
were performed in the mosaic mode, with seven fields covering 
clump 1 in a hexagonal pattern with a central field. The 3 mm re- 
ceivers were tuned at the '^^CO(l-O) frequency (llOGHz), and 
the 1 mm receivers targetted H2CO at 218 GHz. The aiTay con- 
figurations, UV coverage and 1 mm observations are discussed 
in detail in iLeurini et al.l (|2009[) . The receiver temperatures at 
3 mm were around 200 K or better. 

Although continuum emission from clump 1 is not de- 
tected at 1 mm, weak emission is detected at 3mm. This sug- 
gests that the density profile of the clump is rather smooth, 
and that its emission is mostly filtered out by the interferome- 
ter at 1mm. The integrated intensity measured in a 10" diame- 
ter aperture centered on the clump is 0.043 Jy, and 0.12 Jy in a 
20" diameter aperture. Assuming Tdust = 45 K ( see secti on|3.2|l, 
= 2 and /('230GHz = 3.09 x 10 ' cm^ g ' dOssen kopf & H ennind 
'1994), we derive Nh, = 1.6 x 10^^^ cm^^ (resp 1.1 x lO^^ cm^^) 
averaged in a 10" (resp 20") area centered on the clump. 
These values are intermediate and consistent w ithin a factor 
of 2 both with the colu mn densities derived by Leuri ni et al.l 
(2006), and by lLis & S chilke (2003) from H'^CN observations 
(2.6 X 10^^ cm"^). In the following, we will use our newly de- 
rived value Nh, = 1.6 x lO^-' cm"^. 
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3.1.2. Clumps 

Clump 3 was not targetted by the Plateau de Bure observations 
presented in sec. 13.1.11 No continuum was either detected at 
3 mm in the study of iLis & Schilke (2003), which points to the 
fact that the density profile of the clump may be again rathe r 
smooth. From the H'^CN observations of Lis & SchiUcd (l2003h . 
the H2 column density is 1.9 X 10^^ cm"^. 

Using the IRAM 30 m telescope, we targetted the 1-0 and 2- 
1 transitions of C'^O towards this clump. The hyperfine structure 
of the 1-0 transition is clearly resolved. The two observed lines 
have a rotation temperature of 12 K, and lead to a Nh, = 8 x 10^^ 
cm~^, averaged on the extent of the clump (8", see below). 

In the following, we therefore assume an intermediate value 
between the column density as measured from H'^^CN and C'^O, 
= 1.3 X 10^^ cm"^, averaged on the 8" clump. 



the model. Finally, several velocity components, which are sup- 
posed to be non-interacting, can be introduced. 

Table 2. Best fit model results from the CH3OH analysis to- 
wards the two clumps; the uncertainties correspond to the 3cr 
confidence level. The methanol column density (N) and frac- 
tional abundance relative to H2 (x) is computed averaged on the 
extent of the clumps. 



Source 


S 

["] 


[cm-3] 


T 

[K] 


NcH,OH 

[cm-2] 


XcH,OH 


clump 1 
clump 3 


10 

8 


10*^ 
5+^ 10* 


45+47 

35+" 


3+1 10'* 
3+3 10'* 


2x 10-' 
2x10-' 



-5''24'40' 



« -6°25'20' 




R.A. [J2000] 



Fig. 1. In grey scale the H^-'CN (1-0) tr ansition observed wit h 
the Plateau de Bure interferometers by Lis & SchiUcd (l2003l) . 
The numbers indicate the clumps identified by the same authors. 
The solid and dashed circles outline the beams of the APEX 
(~ 290 GHz) and IRAM 30m (~ 241 GHz) telescopes respec- 
tively. 



3.2. Physical conditions in tine clumps : CH3OH analysis 

To determine the properties of the gas in the clumps, we anal- 
ysed the methanol emission at 241 and 290 GHz. The beam of 
the APEX telescope at 290 GHz is almost twice as large as the 
one of the IRAM antenna, and therefore samples different gas 
volumes (see Fig.[Tli. We therefore smoothed the HERA data to 
the resolution of the APEX telescope at 29 GHz. For the anal- 
ysis, we used the technique described by iLeurini et al. I (12004 
for the study of multi-line CH3OH observations, which consists 
of modelling all the lines simultaneously with a synthetic spec- 
trum computed using the Large Velocity Gradient approxima- 
tion, and comparing it to the observations. Rest frequencies are 
from lXu & LovasI ( 1997h. while th e collisional rates were com- 
puted by Pottage e t al.l ( 2002ll20(j4h . The parameters defining the 
synthetic spectrum are: source size, kinetic temperature, column 
density, velocity width and velocity offset (from the systematic 
velocity of the object). The line width and the velocity of the 
object are not free parameters, but are given as input values to 



For our analysis, we modelled the emission towards each 
clump with a single component model, and neglected effects 
due to infrared pumping. We used as free parameters the col- 
umn density of the two symmetry states of methanol, CH3OH-A 
and CH30H-£', the kinetic temperature of the gas, and the H2 
density. For clum p 3, we used a sou rce size of 8" as derived 
from H'^CN by L is & Schilki (|2003|) . For clump 1, the source 
size derived from the analysis of the H'^CN emission is ~ 7"; 
however, other clumps partially fall in the ~ 20" beam size of 
our observations (see Fig. [TJ. Therefore, we adopted a source 
size of 10" which should take into account for the emission of 
the other clumps at the edge of the beam. This corresponds to 
the assumption that all clumps in the Bar have similar physical 
properties. 

The methanol spectra from the two clumps are very similar, 
although clump 3 is shifted in velocity of 0.7 km s-' with respect 
to clump 1 (vlsr ~ 10.0km S-' at clump 1, vlsr ~ 10.7kms-' 
at clump 3). The spectra are characterised by narrow lines 
(~ 1.2kms-'), and no emission is detected in transitions with 
upper level energies corresponding to more than 85 K. As al- 
ready discussed by Leurini et al. (2006), at the original veloc- 
ity resolution of the APEX observations (Av~0.12kms-'), the 
CH3OH ground state lines toward clump 1 show a double- 
peaked profile, probably due to the different clumps sampled 
by the beam. Since no double-peaked profile is detected in the 
IRAM data (Av~0.38kms-') towards the same position, we 
smoothed the APEX spectra towards clump 1 to the same ve- 
locity resolution as the IRAM data. No double-peaked profile 
is detected towards clump 3. However, given the low signal to 
noise ratio of the APEX observations, we smoothed these data to 
a resolution of 0.5 km s-'. The main difference in the methanol 
spectra of the two clumps is in the 1 1 — » Iq-A line at 303.4 GHz 
which is much stronger in clump 3. Since this line has a lower 
level energy of 2.3 K, this suggests that clump 3 is colder than 
clump 1. 

In Table|2]we report the results of the best fit models for both 
clumps. For the column density, the total value of the column 
density for the two symmetry states is given. The high value for 
the spatial density suggests that the gas is close to thermal equi- 
librium, as verified by the results of our LVG calculations that 
show that most lines are partially sub-thermal. We carried out a 
analysis to derive the uncertainties on the kinetic temperature, 
density and column density. The values of the 3cr confidence lev- 
els of each parameters are reported in Table |2] 

All lines are well fitted by the model, with the only exception 
of the 5 1 — > 4i-E transition. The behaviour of this line remains 
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unclear: in massive star forming regions (iLeurini et al.l200'7h the 
5i — » 4i-£' is observed with an intensity roughly half of that of 
the blend of the 5+2 — » 4±2-£ lines, as expected since the three 
transitions have very similar energies. On the other hand, the 
5i — > 4i-E transition in the Orion Bar shows the same intensity 
of the 5±2 — > 4-t2-£ lines blend, and cannot be fitted by our 
models. 

For clump 3, the intensities of the ground state lines in the 
6(1 5it band, and in the li Iq-A line at 303.4 GHz are un- 
derstestimated. The li Iq-A transition is very low in energy, 
and it is expected to be more intense in cold regions. This sug- 
gests that a two component model, with a layer of gas at a lower 
temperature, could be more appropriate for clump 3. However, 
the reduced for our model is 1.8, which is still reasonable for 
a simplified analysis. 

The spectra for clump 1 and clump 3 are very similar. We 
therefore present in Figure|2]the synthetic spectra for the best fit 
model overplotted on the observed spectra for clump 3 only . The 
second feature in Fig. 2(c) is identified with (42,2 — > 32,i) H2CO 
from the lower side band. 



4. Deuterium fractionation in the clumps 

Mapping of the DCN(4-3) emission across the Orion Bar with 
the APEX telescope sho ws that DCN emission originates in the 
H'^CN clumps imaged bv'L is & Schilkel(l2003l) . The overall dis- 
tribution of DCN in the Bar will be the scope of a forthcoming 
paper (Parise et al. in prep). Here, we aim at studying the chem- 
istry at work in the clumps, and therefore we target selected tran- 
sitions of deuterated molecules towards clump 1 and 3. 

In the following subsections we study the DCN excitation 
towards both clumps, present the detection of DCO^ and HDCO 
towards clump 3, and upper limits on the other molecules. 



4. 1 . DCN and H^" CN excitation 

Figure [3] (resp. lU shows the spectra of the several DCN (resp. 
H'-^CN) transitions observed towards the clumps with the IRAM 
30m and the APEX telescopes. In the following subsections we 
study the excitation of these two molecules, using LTE and LVG 
methods. 



4.1 .1 . LTE analysis 

Although the critical densities of the DCN and H'^CN levels 
are quite high, several 10^ to several lO^cm"^, implying that 
the DCN and H'^CN level populations may not be in LTE, it 
is instructive to draw rotational diagrams for the two groups of 
clumps. We corrected the observed line intensities for beam di- 
lution effects by assum ing a size of 7" for clump 1, and 8" for 
clump 3 (as derived bv lLis & Schilkell2003h . The obtained rota- 
tional diagrams are presented in Figure |5] and the derived rota- 
tional temperatures and column densities are listed in Table |4] 

Under the assumption that the lines are optically thin, which 
is confirmed by the hyperfine structure analysis (see section 
I4.I.2I 1. the deviation from linearity in the rotational diagrams can 
be due either to the very different beam sizes of the different 
observations, or to non-LTE effects. The first effect is not com- 
pletely taken out by correcting for beam dilution, because the 
source maybe be actually more extended than the smallest beam. 
This could explain in particular why the DCN(3-2) transition 
(observed in the smallest beam) is weaker than the other transi- 
tions. Although the clumps are very dense (~ 5 x 10^ cm"^), the 



CLUMP 3 

a) 




241.7B8 241.791 241.794 
GHz 



A ., I , ■ , 



ITT 



A 



IT"!!!" '1 



my 



2.4185 10 
Rest Frequency (MHz) 




290.109 290.112 
GHz 



^ J , I ^ . I , I , , '-^^ I , 1 I L^S-l 




2.901 10 2.9015 10 2.902 10 2.9025 10 2.903 10 

Rest Frequency (MHz) 




3.0337 10 3.03375 10 3.0338 

Rest Frequency (MHz) 

Fig. 2. Methanol spectra toward clump 3. a) 5^4 band observed 
with the IRAM 30 m telescope, b) 6— >5 band observed with 
the APEX telescope, c) li-lo-A line (APEX). The synthetic 
spectra resulting from the fit are overlaid in black. The panel 
inserts show the k - Q-A lines in the two bands to illustrate how 
well the models fit the observed spectra. 



critical density of the different levels of the molecules is so high 
that the levels are not populated in LTE. H'^CN seems even fur- 
ther from LTE than DCN, which is consistent with the fact that 
critical densities for H'^CN are around a factor two higher than 
for DCN. The departure from LTE increases with the increas- 
ing upper level energy. The DCN(5-4) transition in clump 1 is 
in particular subthermally excited compared to the three lower 
energy transitions. This line is however blended with an uniden- 
tified line coming from the image sideband (at 350.554 GHz). 
The unidentified line is shifted away from the DCN(5^) lines 
towards clump 3 due to the slight velocity difference between 
the two clumps. 

The square in Figure |5]represents the transition if all flux is 
considered to come from the DCN(5^). This point is still too 
low compared to the three other transitions. It appears that the 
blending of the line cannot account for the low intensity of the 
DCN(5^) emission, and that the J=5 level population is sub- 
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Table 3. Observational results. Fluxes, FWHM and Vki are computed through gaussian fitting. 







clump 1 






clump 3 




Position 




(0", 0") 




(- 


■50", -40") 




Iransition 


J imb tlV 


r W rllVl 


V Lsr 


r T rl^7 

J imb UV 


r W rllVl 


Vlsr 




(Kkms"') 


(kms-^) 


(kms-^) 


(Kkms-^) 


(kms-^) 


(kms-^) 


UV^iNI^Z— i ) 


z.zu±u.u^ 




1 O O-i-O 1 
iU.U±U. 1 


A 1 1 _i_0 O^ 


1 .J i^nis ) 


1 O 7-I-0 1 
lU. /±U. 1 


■nr^vr/a '-)\ 
Lfy^iS V J— Z ) 


1 /^'7_i_0 O^ 


1 ^-i-O 1 
i .0±U. i 


1 O O-I-O 1 
1U.U±U. i 


1 70-I-0 1 1 
1. /U±U. 1 1 


1 1 -i-O 1 


1 O A-i-O 1 
1U.0±U. 1 


DCN(4-3) 


1.15±0.07 


1.3±0.1 


10.1+0.1 


1.75+0.23 


1.4+0.3 


10.8+0.1 


DCN(5^) 


0.25+0.04* 


1.9* 


~ 


0.63+0.10 


1.7+0.1 


10.7+0.1 


H'3CN(l-0) 


1.51±0.05 


1.8 (hfs) 


10.1+0.1 


1.72+0.05 


1.6 (hfs) 


10.6+0.1 




A AQ-i-O OA 


1 1 -t-O 1 


lU.UitU. 1 




1 ^-t-O 1 
1 . j±U. 1 


1 n 7-1-0 1 


H'3CN(4-3) 


1.66±0.06 


1.8±0.1 


10.1+0.1 


1.72+0.10 


2.0+0.1 


10.9+0.1 


DC0+(2-l) 


< 0.03 (3tr) 


1.8* 


- 


0.12+0.01 


1.21+0.14 


10.7+0.1 


H'^^CO+(1-0) 


0.40±0.02 


1.75+0.12 


10.1+0.1 


0.50+0.02 


2.20+0.11 


10.5+0.1 


HCO+(1~0) 


27.8" 


2.08+0.26 


9.9+0.3 


35.3** 


2.17+0.03 


10.5+0.1 


HDCO(2„-lio) 








0.041+0.005 


1.20+0.16 


10.4+0.1 


HfCO(2„-lio) 








0.216+0.014 


1.91+0.15 


10.8+0.1 


C2D(2-1) 








< 0.042 (3a-) 


1.5* 




CH2DOH(2-l) 


< 0.008 (3o-) 


1.8* 




< 0.011 (3a-) 


1.5* 




HDO 








< 0.093 (3a-) 


1.5* 




DNC(2-1) 








< 0.019 (3a-) 


1.5* 




HNC(4-3) 








3.03+0.06 


2.1+0.1 


10.8+0.1 



'Line blended with a line coming from the image sideband. Flux was computed with a two-component gaussian fit, keeping the linewidth fixed. 
** Integrated intensity computed without a Gaussian fit (the line is found to be non-Gaussian). 
*Assumed width to compute the upper limit on the flux. 

Uncertainties given on integrated intensities are just the errors given for a gaussian fit, and do not include the calibration uncertainties (assumed to 
be of the order of 15%). 

All 3a- upper limits are computed using the following relation: J T„,fc dv < 3 x rms x VFWHM x ^v 



Table 4. Results of rotational diagrams and LVG analysis. 



Rotational diagram 

Source T„,(DCN) Ndcn T,.o,(h13cN) Nhlicn DCN/HCN* 

(K) (lO'^cm-^) (K) (10'3cm-2) (%) 

Clump 1 14.0+1.1 1.4+0.3 10.5+0.4 3.1+0.4 0.7+0.2 

Clump 3 11.9+0.5 1.9+0.3 9.9+0.4 2.5+0.3 1.1+0.2 



LVG analysis 



Source 


Tun [DCN] 


Ndcn 


Tk,„[Hi3cN] 


Nhi3cn 


DCN/HCN* 




(K) 


(lO'^cm-^) 


(K) 


(10'-'' cm--) 


(%) 


Clump 1 


18+^ 


1 4+0.4 




8.2!|« 


0.3+0.1 


Clump 3 


32113 


1 0+"'* 


26*4 


1 Q+0.4 
^•'-0.3 


0.8+0.3 



Given error bars are lo". We note that the Ty„ derived from LVG analy- 
sis of DCN and H'^^CN are consistent within the 3cr uncertainty range 
with the Tijii, values derived from methanol. 
•Assuming HCN/H"CN = 70 (Wilson 1 9991) . 



portant for the lower transitions, DCN(2-1) and H'^CN(l-O). 
From the relative intensity of the hyperfine components, it is 

= possible to derive the opacity of the transition. This was done 
using the hfs method of the CLASS software (GILDAS pack- 
age, http://www.iram.fr/IRAMFR/GILDAS/). This method fits 
the spectrum assuming the same excitation temperature for each 
hyperfine component. The results of the fits are shown in Table|5] 
The total opacity of the DCN(2-1) is found to be low. In 
the case where the opacity is well enough constrained so that 
the Tex can be accurately derived, the excitation temperature is 
rather low (Table ID . It should be however noted that the derived 
Tex value is a lower limit because of beam dilution. The excita- 
tion temperature provides a lower limit to the kinetic tempera- 

" ture, as the level populations might not be thermalized, due to 
the high critical densities of each transition. The same holds for 
H'^CN(I-O). The transition is optically thin on both positions, 
and the excitation temperatures are also rather low. 



thermal. After studying the hyperfine structure of the lower rota- 
tional transitions, we will model the line emission using a non- 
LTE method. 

4.1.2. Hyperfine Structure 

DCN and H'^CN rotational transitions have an hyperfine struc- 
ture, caused by the interaction of the electric quadrupole mo- 
ment of the N nucleus (1=1) with the molecular field gradient. 
This causes the transitions to be split into several components, 
reducing the opacity at the line center. This effect is most im- 



4.1.3. DCN and H'^^CN column densities 

In order to derive column densities, we used a standard LVG 
code. The hyperfine structure of the lower levels is not ex- 
plicitely taken into account. However, to account for the re- 
duced opacity due to the hyperfine splitting, we replaced the 
escape probability by a weighted mean of the escape proba- 

v-i 1 - e-f''"'"' 

bili ties of the hyperfine components: /3 - / fi .For 

^ fiTtot 

both isotopolog ues, we used co l lision coefficients calculated 
for HCN-He by IWernli & Faurg (|2009|) . These collision coef- 
ficients were uniformly scaled by a factor 1.37 (square root 
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Table 5. Results of the hyperfine structure fit. 



Source 


Line 




FWHM 


Tmb '''tot 




T 

^ ex 


n,,-(50K) 






(km s-') 


(km s-') 


(K) 




(K) 


(cm-3) 


Clump 1 


DCN(2-1) 


10.01±0.02 


1.90±0.04 


1.2±0.1 


0.7±0.2 


5.0±1.8 


9.1x10^ 




H"CN(l--0) 


10.13±0.02 


1.78±0.07 


0.92±0.09 


0.56±0.41 


4.6+3.8 


1.8x10^ 


Clump 3 


DCN(2~1) 


10.66±0.01 


1.30±0.02 


2.8±0.1 


0.1±0.4 


31.5±127 


9.1x10^ 




H'3CN(l-0) 


10.62±0.02 


1.61±0.03 


0.96±0.02 


0.10±0.04 


12.7±5.3 


1.8x10^ 



Clump 1 f 
(0",0") position / 


\ 

^ DCN(3-2) 





Clump 3 
(-50", -40") position 



A DCN(3-2) - 



DCN(8-1) 



DCN(4-3) 




i tcN(5-4)x2 



-10 10 20 30 

Velocity (km/s) 

Fig. 3. DCN spectra observed towards the two groups of clumps. 
Hyperfine structure fits are displayed in green for the DCN(2-1) 
lines. A parasite line coming from the image sideband is visible 
in the DCN(5^) data (clump 3). Due to the slight velocity dif- 
ference between the two clumps, the parasite line is blending the 
DCN(5-4) fine on clump 1. 



of the ratio of the reduced mass of the HCN-He and HCN-H2 
systems) to approximate HCN-H2 collisions . This approxima- 
tion is generally supposed by astronomers to be valid for HCN- 
pH2 collisions, as para-Ha has a spherical symmetry in its J=0 
state. Theoretical computations have shown that the accuracy 
of this assumption is however limited, the discrepancy for in- 
dividual rates reachin g up to a fa ctor of two for p articular sys- 
tems (H^O-H./H cPhillips et al.lll996.. CO-H^/He.lWernh et alJ 
20061: HC^N-H9/He. lWemli et al.l2007l: SiS-H7/H e. Lique et al. 



2008 1). or even higher (NH^-H7/He. iMaret et al.l l2009)l This is 
however the only solution as long as HCN-H2 collision rates are 
not available. 

Input parameters for an LVG code are nn,, T and Nmoi/Av. 
Comparing modelled line temperatures with observations brings 
another parameter into play, the beam filling factor. On the ob- 
servational side, we want to fit simultaneously for each molecule 
the observations of all available transitions and the opacity de- 
rived by the hfs method for the lower transition. There are thus 



, e -2. - 



Clump 1 J 
(0",0") position [ 



1 lU H^^CN(l-O) x5 




Clump 3 
(-50", -40") position 



^1 r^j^R — 




10 20 

Velocity (km/s) 



Fig. 4. H'^CN spectra observed towards the two groups of 
clumps. Hyperfine structure fits are displayed in green for the 
Hi^CN(l-O) lines. 



5 observables for DCN and 4 for H'^CN. It is statistically rather 
meaningless to try to fit a sample of four/five observations with a 
4-parameter model. We therefo re use i ndependently determined 
information for the source size dLis & Schilke 2003), and the H2 
density as derived from methanol observations (see section 3.2), 
and only fit the temperature and the DCN and H'^'CN column 
densities. 

Figure |6] shows as an example the x~ analysis for fitting 
LVG models for the DCN molecule towards clump 1 . The 
was computed using the 4 DCN transitions, weighted by uncer- 
tainties including a 15% calibration uncertainty of the integrated 
flux, as well as the opacity of the DCN(2-1) line as derived in 
14.1.21 The derived kinetic temperature is consistent within the 
3cr uncertainty range with the one derived from the methanol 
analysis. The DCN column density compares also very well with 
the rotational diagram result. The Icr confidence interval for the 
two parameters is shown in Table |4] We analyzed in the same 
way the DCN emission towards clump 3 as well as the H'^CN 
towards both clumps. The results are all shown in Table |4] 
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Fig. 5. Rotational diagrams for DCN (triangles) and H'-'CN (dia- 
monds). The square represents the upper limit for the DCN(5^) 
integrated intensity towards clump 1, computed as the total inte- 
grated intensity including the parasite line (see text). 
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Fig. 6. LVG analysis results for the DCN molecule towards 
Clump 1. 



4.2. Low DCO"^ emission 

DCO+(2-l), H'^CO+(1-0) and HCO+(1-0) were observed to- 
wards the two clumps. The spectra are presented in Figure|2] By 
comparison with H'^CO+(1-0), HCO+(1-0) is found to be op- 
tically thin. The observed parameters are listed in Table [3] Note 
that toward clump 3 the linewidth of the H'^CO+ and HCO+ 
lines are significantly larger than the HCN and DCN linewidths 
(as determined from the hyperfine structure fit). The thermal 
broadening expected for H'^CO^ at 50K is 0.27 kms ', and 
0.47 km s"' at 150 K, the temperature characteristic of the in- 
terclump gas. This points to the fact that HCO^ and H'^CO^ 
are also present in the warmer interclump gas, in agreement 
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(0",0") position DC0^(2-1) 
Lvv>ir-n-'-''V-''\juJui/^^ 




HCO (l-0)/70 



Clump 3 

(-50", -40") position 




10 20 
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Fig. 7. DCO+, H'^CO^ and HCO+ spectra. 



with the interferomet ric studies of I Young Owl et all (l2000h and 
& Schilkd (I2003h . who claimed that a large fraction of the 
emission of H'^CO^(l-O) is extended and thus filtered out by 
the interferometer. 

On the other hand, we do not expect deuterated molecules to 
be at all present in the hot interclump gas in steady-state. The 
fact that the main isotopomers are likely to trace both the two 
gas components implies that the D/H ratio derived from the ob- 
servations should be considered as a lower limit. 

We estimate column densities towards the clumps (assuming 
resp. sizes of 7" and 8" for clump 1 and 3) for diff'erent rota- 
tional temperatures, assuming that all emission comes from the 
clumps. Although DCO^ and H'^CO^ column densities vary by 
a factor of two to three for temperatures between 20 and 70 K, 
the abundance ratio is less affected by the temperature uncer- 
tainty (less than a factor of 1.5, see Table |6]l 

Can we quantify how mu ch of the HCO^ e mission actu- 
ally comes from the clumps? iHogerheijde et al.l (|l995) mod- 
elled multi-frequency transitions of several molecules with a 
two-phase (interclump + clump) model, and they concluded that 
10% of the column density of the material is in the clumps. We 
can then get an estimate of the clump contribution to H'^CO^ 
emission. Since the emission is likely dominated by the inter- 
clump gas, we estimate the column density of HCO^ assuming a 
temperature of 150 K. For clump 3, we then find a column den- 
sity of H'''CO+ of 6.2x10'^ cnr^. If the clump represents 10% 
of the column density, then N(H'^CO+) ~ 6x10'^ cm - in the 
clump, which is roughly a factor of two to three lower than what 
we estimated. We may thus underestimate the DCO^/ HCO^ by 
a factor up to three. 
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clump 1 






clump 3 






Trot 


Noco+ 


Nh13co+ 


Ndco+/Nhco+' 


Ndco+ 


Nh13co+ 


Ndco+/Nhco+ 


1 




(cm-2) 


(cm"-) 




(cm"-) 


(cm"-) 






20 K 


<1.6xl0" 


1.0x10'^* 


< 2.3x10"" 


(5.0±0.8)xlO" 


(9.9±0.1)xl0'2 


(7.3±1.4)xl0- 


-4 


30 K 


<2.0xl0" 


1.4x10'^ 


< 2.0x10"'* 


(6.2±1.0)xl0" 


(1.4±0.2)xl0" 


(6.4±1.3)xl0- 


-4 


40 K 


<2.4xl0" 


LSxlO'"* 


< 1.9x10""* 


(7.5±1.2)xl0" 


(1.8±0.2)xl0'3 


(6.0±l.l)xlO- 


-4 


50 K 


<2.8xl0" 


2.2x10'^ 


< 1.8x10"-* 


(8.9±1.4)xl0" 


(2.2±0.2)xl0" 


(5.9±1.0)xl0- 


-4 


70 K 


<3.7xlO" 


3.1x10" 


< 1.7x10"-* 


(1.2±0.2)xl0'2 


(3.0±0.3)xlO" 


(5.7±l.l)xl0- 


-4 



'We assume HCO+/H'''CO+ = 70 jWilsodl 19991) . 



4.3. Detection of HDCO 

In order to test the possibility that the deuteration in the clumps 
originates in chemistry involving CHaD^, we targetted HDCO, 
a molecule synthesized in the gas phase v i a the CHiD"^ route, 
through the following reactions dTurneJ 120011 : iRoueff et alj 
l2007h : 

CH2D+ + H2 ^ CH4D+ (4) 
CH4D+ + e^ CH2D + H2 (5) 
CH2D + O ^ HDCO + H (6). 



derived HDCO/H2CO ratio above is thus to be considered as a 
lower limit (see discussion in section |53] l. 



4.4. Upper limits on other interesting deuterated moiecuies 

In the following, we derive upper limits on the fractional abun- 
dances of other deuterated molecules that were not detected. The 
discussion of the results is left to section|5] 



4.4.1. C2D 
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Fig. 8. HDCO and H^^CO detected towards Clump 3. 



C2D is another molecule believed to form in the gas phase 
from CH7D+. RoueflF et al. (2007) predicted an D/H ratio for this 
molecule of 3.9x10"^, and an abundance of 6.6x10"" with re- 
spect to H2, in their low metal model at 50 K. 

The J = 2-1 band was observed towards clump 3. The 
VESPA correlator at the IRAM telescope showed unfortunately 
a lot of platforming. Using the IMHz filter bank, the lines are 
not detected, and the rms is 8mK {T,„b scale) at a resolution 
of 2.1 kms"'. Assuming a linewidth of 1.5 kms"' (as found in 
HCN and DCN observations), we derive the flux upper limit 
listed in Table|3] Assuming a Tiot of 35 K, this corresponds to an 
upper limit on the C2D column density of 2.5 x 10'^ cm"^, and 
an upper limit for the fraction al abundance of 2 x 10"'", compat- 
ible with the the prediction of IRoueff' et"aLl (|2007|) . 



4.4.2. DNC 



We detected the HDCO(2-l) transition towards clump 3 
(Fig.O. The line has a Vki of 10.4 kms"', consistent with other 
lines detected towards this clump. The line width (1.2 kms"') is 
also similar to those of other deuterated species (DCN, DCO^), 
making the detection of HDCO rather secure. We also detected 
the H2'^CO(2-l) line towards the same position. 

We computed the column density of the two molecules us- 
ing the LTE approximation. The column densities vary by a fac- 
tor of three depending on the assumed temperature between 20 
and 70 K. However, as the two observed lines have similar en- 
ergies, the HDCO/H2'^CO ratio is less sensitive to the tempera- 
ture, varying only by a factor of 1.2 over the whole temperature 
range, which lies within the error bar of the determined ratio. We 
give in Table [T] the column densities and D/H ratio for a fixed 
temperature of 35 K (as pointed by the analysis of the CH3OH 
emission). We assume again H2'^CO/H2CO = 70 (Wilson 1999). 

One shou ld note that H2CO is also partly tracing the inter- 
clump gas (iLeurini et al.l l2006l Leurini et al. submitted). The 



DNC is synthesized in the gas phase, mainly from a route involv- 
ing the H2D^ ion, as opposed to DCN which can be synthesized 
from CH2D+ ( Turner 2001: Rou es' et al. 2007.). We searched for 
the 2-1 transition of this species, to get constraints on the im- 
portance of the H2D^ chemistry in clump 3. We did not de- 
tect it, at a rms noise level of 9 mK (T,,,/, scale) and a resolu- 
tion of 0.3 1 km s"' . We derive the flux upper limit listed in Table 
|3] This corresponds to an upper limit of 1.5x10" cm"^ for the 
DNC column density (assuming Tim = 35 K), and an upper limit 
on the abundance relative to H2 of 1.0x10"'^. Assuming that 
the HNC(4-3) line is optically thin, we find Nhnc = 1.05x10'^ 
cm^. This has to be considered as a lower limit if HNC(4- 
3) is not optically thin. This translates into an upper limit of 
DNC/HNC < 1.4x10"^ (3cr). This upper hmit is not very con- 
straining compared to the deuteration ratios measured in HCN, 
because HNC is found to be more than three orders of magnitude 
less abundant than HCN. The very high HCN/HNC ratio in the 
clump is to our knowledge one of the highest observed so far. 
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4.4.3. CH2DOH 

We searched for monodeuterated methanol, which is believed 
to trace the evaporation of highly deuterated ices ('Parise et al.' 
[2002, 2004, 2006). We looked for the 2k- Ik rotational band 
at 89 GHz, corresponding to low energy transitions. CH2DOH 
was not detected towards any of the two clumps. Rms levels of 
4 mK toward clump 1 and 6 mK toward clump 3 (T,,,;, scale) were 
reached, at a 0.26 km/s resolution. Assuming a linewidth of 1.8 
(resp. 1.5) kms"' for clump 1 (resp. 3), we derive the upper limit 
for the integrated intensity listed in Table[3] This corresponds to 
an upper limit of 1.9 x 10'"* cm"^ on the CH2DOH column den- 
sity in clump 3, i.e. an upper limit on the fractional abundance 
of l.SxlO-''. 

4.4.4. HDO 

The HDO line at 241 GHz was searched for towards clump 3 but 
not detected, at a noise level of 47 mK (rms, T,„h scale), and a 
resolution of 0.3 kms The associated upper limit for the in- 
tegrated line intensity given in Table [3] corresponds to an upper 
limit of 4.4 x 10'^ cm"^ for the HDO column density, i.e. an up- 
per limit for the abundance of 3.4x10"'". 



5. Discussion 

The main results of this study are the confirmation through 
multi-transition observation s of the significant deuterium frac- 
tionation of HCN found by iLeurini et al.l (l2006l ). as well as the 
detection of significant fractionation of HDCO in dense clumps 
in the Orion Bar PDR that are too warm for deuteration to be sus- 
tained by the H2D^ precursor. In the following, we consider the 
possible explanations for this high deuteration: products evapo- 
rated from ices surrounding dust grains, or gas-phase products. 
This study can also help understand the chemistry at work in 
PDRs. 

5.1. Grain evaporation ? 

One plausible origin for highly deuterated molecules is the 
evaporation of ices surrounding dust grains. Highly deuter- 
ated methanol, formaldehyde and water are observed in hot 
core s (e.g. [jac q et al. 1988; Comito et al. 2003J and hot cori- 
nos (IParise et al.(l2002l |2004 l2005l l2006l) . and are thought to be 
remnants of the cold prestellar phase. 

In particular methanol is believed to form primarily on dust 
grains, as the gas-phase production routes for this molecule are 
not efficient. The same may also be true for water, the main 
constituent of ices surrounding dust grains. Deuterated isotopo- 
logues of these molecules may thus trace surface chemistry pro- 
cesses that enhance the deuterium content of these species, be- 
cause the accreting atomic D/H ratio is much higher than the 
elemental D/H (~ 10"^). Atomic accretion on dust grains and 
surface reactions take place in cold environment, i.e. in condi- 
tions which are also favoring a large atomic D/H. Although it 
is not clear which thermal history the gas in the Orion Bar has 
undergone in the past (are the clumps remnants of high density 
clumps of the cold molecular cloud, which are still shielded from 
the PDR radiation?), it may be possible that we are witnessing 
in the clumps evaporation of ices that were formed during the 
cold molecular cloud phase. In this case, we would expect the 
composition of the ices to be comparable to the ones evaporated 
in hot corinos, i.e. with high CH2DOH/CH3OH ratios. 



The temperature of the grains should follow closely the gas 
temperature (typically <50 K) in the clumps, due to the high den- 
sity in the clumps ( Kruegel & Walmslev 1984). Although this 
temperature is high enough to sublimate CO ices (Tevap~20 K), 
it is not high enough to sublimate polar ices dominated by 
H2O (Tevap~100K). Pure HCN ices are expected to sublimate in 
the ISM conditions at temperatures ~ 65 K (extrapolation from 
the sublimation temperatures in comets from Prialnik 2006). 
However, HCN is only a minor constituent of ices in objects 
where it was dete cted (< 3% of water in W33A, .25% of wa- 
ter in Hale Bopp. lEhrenfreund & Charnlevll2000l) . It will thus 
behave as small impurity in the water ice, and its sublimation 
will be dominated by the sublimation of H2O. It is thus rather 
unlikely that evaporation plays an important role in the warm 
clumps, except maybe on the clump surfaces. 

The low relative abundance of metha nol (2x1 0~^) com pared 
to the abundance in hot corinos (~ 10~^. Maret et al! 2005) may 
be a further sign that thermal evaporation is not the dominating 
chemical process, but rather that gas-phase chemistry is the most 
significant mechanism responsible for the enhanced deuteration 
levels. Note that grain chemistry models including non-thermal 
desorption of methanol are e nough to explain t hese relatively 
low abundances of methanol dGarrod et alj|2007 ). Of course it 
is also possible that evaporated ices get photodissociated in the 
PDR, and this effect will be studied in more detail in a forthcom- 
ing paper (Parise et al. in prep). We focus here now on an alter- 
native explanation, assuming that gas phase reactions are more 
predominant. 

5.2. Gas-pliase cliemistry ? 

The low temperature of the clumps with respect to the subli- 
mation temperature of water ices points to the fact that thermal 
ice evaporation does not play a predominant role in the clumps. 
At the same time, the grain temperature (which, at these high 
densities, should follow the observed gas temperature) is too 
high for CO to stick efficiently on dust grains, preventing ef- 
ficient grain chemistry to take place. The present chemistry is 
thus likely dominated by gas-phase processes, and the clumps 
may be in this sense remnants of the molecular cloud, which 
have been warmed up to temperatures higher than for typical 
molecular clouds, due to their location behind the photoionis- 
ing front. The CH3OH abundance, in particular, is similar to that 
measured towards the TMC-1 dark cloud (3x10"'^, S mith et al.l 
i2004l) . and this molecule appears to be a remnant of grain chem- 
istry happening during the colder er a of the clou d, followed 
by non-thermal desorbing processes (iGarrod et al.l 12007). This 
possibility may be further tested by observing the D/H ratio in 
methanol. A formation of methanol during the cold molecular 
cloud phase is expected to lead to high CH2DOH/CH3OH ra- 
tios, as measured e.g . in prestellar cores (between 5% and 30% 
Bac mann et al. I I2OO7I) . Unfortunately, the upper limit we derive 
here (CH2DOH/CH3OH < 0.8) is not constraining enough, and 
further insight will require significantly deeper integrations, for 
which the sensitivity of ALMA will be needed. 

For other molecules that can also form in the gas-phase, gas- 
phase processes are likely to be predominant. To test the hypoth- 
esis that gas-phase chemistry is the dominating process in the 
clumps, we compare in the remaining discussion the observed 
abundances and D/H ratios to the results of a pure steady-state 
gas-phase model. In particular, we are interested to see if high 
DCN/HCN ratio, low DCO+/HCO+ ratio and the detection of 
HDCO can be explained in the framework of a simple steady- 
state model. 
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Table 7. Summary of column densities, abundances and D/H ratios in the observed molecules (assuming Trot=35 K for species 
detected with only one transition). 



Molecule 


N 

(cm-2) 


clump 1 

X 

(cm-3) 


XD/XH 


N 

(cm-2) 


clump 3 

X 

(cm-3) 


XD/XH 


H"CN 
DCN 


(3.1±0.4)xl0" 
(1.4±0.3)xl0" 


1.9x10-'" 
8.8x10-" 


0.7+0.2 % 


(2.5±0.3)xl0'3 
(1.9±0.3)xl0'3 


1.9x10-'" 
1.5x10-'" 


1.1 ±0.2 % 


H"CO+ 
DCO+ 


(2.0±1.0)xl0'' 
<2.2xl0" 


1.3x10-'" 
< 1.4x10-'^ 


<2xl0-'' 


(1.6±0.2)xl0'' 
(6.9±l.l)xl0" 


1.2x10-'" 
5.3xlO-'2 


(6.1±1.1) 10-* 


HDCO 








(1.2±0.1)xl0''' 
(4.8±0.8)xl0'2 


9.2x10-" 
3.7x10-" 


0.6±0.1 % 


CzD 








<2.5xl0'3 


<2xl0-'" 




HNC 
DNC 








1.1x10" 
< 1.5x10" 


< lxl0-'2 


< 1.4 % 


CHjDOH 


< 1.7x10"* 


< 1.1x10-' 




< 1.9x10'-* 


< 1.5x10-' 




HDO 








<4.4xl0'3 


< 3.4x10-'" 





5.3. Comparison to gas-phase model 

We compare in the following the observed fractionation and rel- 
ative abundances of molecules to the predictions of a pure gas- 
phase model, based on an update d version of the chemical model 
described bv lRoueff et al.l (l2007l) . The new model takes into ac- 
count the new branching ratios of the NaH^ dissociative recom- 
bination (Moleket al. 2007), computes the pre-exponential fac- 
tors of the reverse reactions involved in the deuterium fractiona- 
tion of CH3 with the proper factors involved in the translational 
partition functions. In addition we have included the branching 
ratios of the electronic recombination of HCO^ from Amano 
(I199 O) where the channel towards CO is found to be predom- 
inant. The radiative association reactions of CH^ and deuterated 
substit ues with H2 have been derived from the theoretical predic- 
tions of lBacchus-Montabonel et all (l2000h who display values at 
different temperatures. 

The exothermicity of the CH^ + HD reaction and subsequent 
deuteration steps are not well con strained. These exo thermicities 
were derived experimentally by ISmith et aTl ( Il982t) . However, 
theoretical assessments from zero-point energies lead to higher 
barriers. This has the effect to allow deuteration to remain ef- 
ficient at even higher temperatures (up to 70 K). In the follow- 
ing, we will compare our observations to two models, calculated 
resp. by using the experimental values from ISmith et al.l (1 19821 
sua model) and the exothermicities computed theoretically from 
zero-point energies (theo model). The respective exothermicities 
assumed in the models are listed in Tab. [8] 

Figures ?? and |9] show the observed ratios towards clump 
1 and 3 respectively, as well as the predicted ratios from the 
chemical models. The solid curves correspond to the saa model, 
and the dashed curves to the theo model. The two curves for 
each model correspond to the densities 3x10^ and 10^ cm--^ 
(in order of increasing D/H ratio). The grey filling delimitates 
the observed values (Icr) or upper limits (3cr), and the tem- 
perature range (3cr) as derived from the CH3OH analysis. We 
computed the model predictio ns for the two diff erent sets of 
elemental abundances used in iRoueff et al.l (l2007h . The "warm 
core" elemental abundances are representative of a mostly un- 
depleted gas, while the "low metal" case, shown to lead to the 
best agreement to observations towards dense molecular clouds 



Table 8. Exothermicities used in the two different chemical mod- 
els. 



Reaction 


Exothermicity 


Exothermicity 




saa model 


theo model 


CH+ + HD «^ CH,D+ + H2 


370 K 


670 K 


CH;D+ + HD CHD+ + H2 


369 K 


433 K 


CHD^ + HD CD+ + H2 


379 K 


443 K 


CH+ +D2 ^ CHD+ + H2 


713 K 


1005 K 


CHJ + D2 « CH2D+ + HD 


319 K 


592 K 


CH2D+ + D2^ CD+ + H2 


599 K 


564 K 


CH2D+ D2 « CHb^ + HD 


317 K 


354 K 


CHDJ -H D, CD+ + HD 


290 K 


151 K 



dGraedel et al. !'l982'). involves moderate depletions of C, N and 
O, and strong depletions of S and Fe. 

The models all show a decreasing fractionation with increas- 
ing temperature. In the case of clump 1 (Figure ??), DCO^ 
was not detected, and thus we have only an upper limit for 
the observed DCO+/HCO+ ratio. The observed DCN/HCN ra- 
tio is consistent with the model prediction for a clump tem- 
perature of ~ 45 K, for the saa model with low-metallicity el- 
emental abundances. The theo model tends to overestimate the 
DCN/HCN ratio, requiring temperatures >70K to reproduce 
the observed ratio, which are still consistent but rather at the 
high end of the uncertainty range of the temperature as de- 
termined from the methanol observations. Both models with 
warm core elemental abundances also qualitatively reproduce 
the DCN/HCN for a reasonable temperature range. The upper 
limit on the DCO^/HCO^ ratio is in the case of this clump not 
very constraining. Depending on the model, it points to T > 33 K 
or T > 40 K, which is consistent with the temperatures given by 
the DCN/HCN ratios. 

In the case of clump 3 (Figure IS, both DCO+/HCO+ and 
HDCO/H2CO ratios have been observed. The derived fractiona- 
tion ratios have been estimated assuming that the main isotopo- 
logue is only tracing the clump. However, there is some evidence 
that H2CO and HCO^ are also present in the interclump medium. 
The case for HCO^ was discussed in section |4~2l where we ar- 
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gued that we may have underestimated the DCO^/HCO^ ratio 
by up to a factor of three. H2CO has been shown to be roughly 
as abundant in the clumps as in the interclump gas (Leurini et al. 
submitted). They derive an H2CO column density of 6x10'^ to 
6x10'^ cm""* in the interclump gas, and of (2.8-5.4) 10'"* cm""* 
towards clump 1 , which is up to a factor of three lower than our 
estimate from Hj'^CO towards clump 3. We may thus be overes- 
timating the H2CO abundance in clump 3 by a factor of two to 
three. As a consequence, the ratios we have derived are in fact 
only to be considered as lower limits for the true ratios in the 
clumps (represented by the ascending arrows on Fig. |9]l. Note 
that this is not the case for the DCN/HCN ratio, because H'^CN 
has been shown to trace only the clumps (Lis & SchilkQ l2003h . 
In the case of clump 3, the DCN/HCN ratio points to a temper- 
ature of ~45 K (low-metal elemental abundances), or 20-50 K 
(warm core conditions). Both DCO+/HCO+ and HDCO/H2CO 
lower limits are in principle consistent with the model. However, 
in the case of the low-metal saa model, for a temperature of 
45 K, agreement with the model would require that we increase 
the H2CO abundance in the clump by a factor of six, which 
seems not consistent with the study of Leurini et al. (submit- 
ted). For this clump, better qualitative agreement is found with 
the models with warm core elemental conditions. It is also in- 
structive to see if the upper limits on other deuterated molecules 
(C2D, DNC, HDO) are consistent with the models (Figure [TOl). 
In the case of HDO and C2D, where we have no observations 
of the main isotopologue, we compare the abundances directly 
with the model predictions. C2D upper limits are not constrain- 
ing — a two orders of magnitude increase in the sensitivity would 
be required to add interesting constraints. The DNC/HNC upper 
limits point to T > 30 K. But the most constraining molecule is 
certainly HDO, which non-detection is consistent with the mod- 
els only for T> 40-50 K. The fact that the HDO abundance seems 
to be consistent and even rather on the lower side of the model 
predictions is certainly a good hint that thermal evaporation is 
not playing a dominant role in injecting water molecules into the 
gas-phase. 

As a conclusion from this comparison, we find a good quali- 
tative agreement between the observations and pure steady-state 
gas-phase models. Although such a comparison is obviously 
limited, because the physical structure (temperature and density 
gradient) of the clump is not taken into account, and because 
the chemistry is only computed at steady-state, it shows that the 
scenario of warm gas phase deuterium chemistry is viable for 
explaining the high deuteration ratios observed in DCN towards 
the Orion Bar clumps. 

The comparison with the models is unfortunately limited 
by our poor handle on the clump temperature. Good probes 
of the temperature are usually the inversion l ines of the NH3 
molecule. NH3 was mapped in the Orion Bar bvl Batrla & WilsonI 
( 1200 3) with single-dish. Although they find that the NH3 emis- 
sion comes mostly from the same regions as the HCN emission, 
they find high kinetic temperatures (T> lOOK) and argue that 
ammonia is located in the surface layers of the clumps, where 
icy mantle around dust grains evaporate. The traced temperature 
is thus rather typical of that of the interclump gas, and NH3 is 
therefore not a good probe of the clump temperature. Deeper 
observations of several lines of CH3OH, including high- J transi- 
tions might allow to constrain the temperature better. 

On the model side, a more detail treatment taking into ac- 
count photodissociation in the PDR is the next step to take. 



6. Conclusions 

We have presented observations of deuterated molecules towards 
two dense clumps in the Orion Bar photodissociation region. 
These observations were aimed at confirming and understand- 
ing the origin of the DCN emission first detected in this region 
bv Leurini et al.l (l2()06i) . We confirmed through the observation 
of four transitions of DCN the detection of this molecule towards 
one clump, and detected it towards a second clump of the Orion 
Bar We also detected DCO^ and HDCO towards this second 
clump, and provided upper limits for the abundance of other rel- 
evant deuterated molecules. From the observation of these sev- 
eral species, formed by chemistry induced either by H2D^ or 
CH2D^, we find evidence based on a pure gas-phase chemistry 
model that the main ion responsible for deuteron transfer in the 
Orion Bar is CH2D^, as opposed to previously observed cases 
of colder regions or hot cores where H2D^ was the main ac- 
tor (in the case of hot cores, deuterium fractionation is believed 
to be a fossil of cold chemistry in the earlier cold evolutionary 
phase, preserved into ice mantles). The luke-warm conditions 
in the Orion Bar clumps thus allowed us to observationally test 
chemical models in a different temperature range than most pre- 
vious studies dealing with deuterium fractionation. More refined 
understanding of the chemistry at work in the Orion Bar will re- 
quke more detailed chemical modelling, coupling PDR models 
with gas-grain chemical networks. 
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Fig. 10. Undetected molecules towards clump 3. Same plotting codes as for Fig.|9l 



